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Abstract 

We present some old and new results in the enumeration of random walks in 
one dimension, mostly developed in works of enumerative combinatorics. The 
relation between the trace of the n-th power of a tridiagonal matrix and the 
enumeration of weighted paths of n steps allows an easier combinatorial enu- 
meration of paths. It also seems promising for the theory of tridiagonal random 
matrices . 



1 Introduction. 

Already at the foundation of the theory of random matrices by E.Wigner, the 
relevance of the combinatorics of random walks was recognized {jj . The follow- 
ing decades witnessed an explosion of the theory and the applications of random 
matrices. A number of specific techniques were devised and the relation with 
the combinatorics of random walks was almost forgotten, with few remarkable 
exceptions H. An interesting and recently investigated open problem, ||, Q, 
[fH i O ' @ ' H ' * s ^ ne spectral density of non-hermitian tridiagonal random ma- 
trices. Here the enumeration of one dimensional random walks, where to each 
step an independent random variable is associated, plays a dominant role [jllj . 
Since the subject of random walks is very basic and useful, it has been thor- 
oughly studied both by mathematicians and physicists for a long time, yet with 
techniques and aims so different that the relevant literature about inhomoge- 
neous walks of one group ignores that of the other group. The aim of this paper 
is to review some old and recent results in enumerative combinatorics of random 
walks pertinent to the analysis of tridiagonal matrices in a language accessible 
to physicists, with examples, comments and relations that are not easily avail- 
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able in the literature. 

A simple walk of p steps on the lattice of integers, may be coded by a 
sequence {/ii, ..,/i p }, where each Hk may take the values 1 (right move) or — 1 
(left move). The walk is simple by the fact that we exclude the move fik = 0. 
The index k can be regarded as discrete time. Since our walks start from the 
origin, at time k the walker is on site s(k) = ji\ + .. + fik- It is very convenient 
to view the walk as a continuous broken line in the lattice {k, s) of time and site 
occupancy. Most of the terminology is based on this picture. An illustration of 
a simple walk we consider is given in Fig.l; the discrete time that numbers the 
steps is measured in the horizontal axis, and the positions of the walker on the 
line are recorded in the vertical axis. 




Fig.l 

A little terminology is useful. Simple walks where also the end site is the 
origin, that is s(p) = 0, are "closed" . A simple closed walk is "weakly positive" 
if s(k) > for all k. If we further require that it returns to the origin only 
at the end, the simple walk is " strictly positive" . Weakly positive closed walks 
are also known as Dyck paths. The width of a walk is the number of different 
visited sites. 

Simple walks can be classified and then counted according to different prop- 
erties. A large number of results in the enumeration of simple walks have been 
obtained by a method which we call "iteration by length" : the counting of walks 
with a given property is related to the counting of similar but shorter ones, at- 
tached sequentially. In terms of generating functions, the method leads to a 
system of equations which, in several important cases, can be solved explicitly. 
One of the best known examples is the enumeration of returns to the origin 
in closed positive walks. In sect. 2 we illustrate the method by evaluating the 
number of peaks and valleys in several ensembles of walks. It is a very simple 
exercise given the known methods in enumerative combinatorics [|l2] , |l5| ] , 
[ pT[ , and the method is familiar to physicists as it recalls a Dyson equation in 
quantum field theory. 
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Recently we developed a technique , to count the number of visits of 
each site, for an ensemble of simple walks. The method is again recursive but 
here properties of longer and wider walks are computed in terms of shorter and 
narrower ones. We call this method "iteration by top insertions". It leads to 
explicit results in the form of sums of products of binomial coefficients, the 
sum being actually in several indices taking values in the compositions of an 
integer, occasionally with further restrictions. Since the results are involved, 
techniques to simplify them would be very useful. We think that this method is 
promising in the analysis of the spectrum of tridiagonal random matrices and 
the related random hopping problem in one dimension, where few exact results 
are available, especially in the non-hermitian case. In this introduction we shall 
describe in detail the connection of this counting problem for simple walks with 
tridiagonal matrices. 

We then found that our method 'iteration by top insertions' had already been 



discovered in the form of enumeration of vertices of planar rooted trees , 13 , 
[ pl| , p^j . We present this relation in sect. 3 . 

Occasionally a problem may be analyzed both by the simpler method of 
sect. 2 as well as by the more complex one of sect. 3. We provide an example 
by evaluating the area under Dyck paths, recently discussed by Jonsson and 
Wheater J2j| with the simple method. While our discussion does not add to 
their solution, the interpretation in terms of a tridiagonal matrix facilitates its 
determination. 

The method in sect. 2 is most useful for homogeneous (or almost homogeneous) 
random walks, where all steps (or almost all of them) have the same probability. 
Some early papers are fl5| , |l6| . We remark that even the completely inhomoge- 
neous walks where each step is associated with an independent random variable 
were enumerated [ p"7| by methods similar to the one we describe in sect. 2. 

As it is often the case, computing the same quantities in two different ways 
originates identities which could be hard to prove in a direct way. For example, 
the sums over compositions of an integer of products of binomials introduced in 
sect. 3, will be shown to yield simple trigonometric sums. An early example of 
this type of sums is provided by eq.(22) in a paper by D.A.Klarner fl8f . 

Most of the methods described in this paper have been developed for walks 
where also the move fik = 0, that is an horizontal step , is allowed. Dyck paths 
properly generalized to include such moves are called Motzkin paths. Enumer- 
ation formulas become more complicated and, for the introductory purpose of 
this paper, we avoid these walks. Some results for generating functions for Dyck 
and Motzkin paths and their relation to continued fractions were recently sum- 
marized by C.Krattenthaler Jl9| . 

Ensembles of walks are among the most basic topics in probability theory 
and their statistical properties are relevant in a large class of models in statis- 
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tical mechanics. 

In one dimension, to each step eu of a walk 7 one may assign a weight w(ek)- 
The whole path is then attributed a weight 10(7) = w{e{]w{e2)--w{e p ). Finally 
one is interested in summing over walks in some set T, weighted with their own 
factor w(j). One encounters this procedure in the discrete approximation of 
functional integrals J T>x f[x(t)], where the integration over continuous positive 
functions x(t) with boundary conditions x(t ) = = x(tf) is replaced by an 
ensemble of Dyck paths. The enumeration methods described here may be use- 
ful for the evaluation of functional integrals where the paths may have various 
restrictions. In the example 3 of sect. 3, few cases of weighted sums of Dyck 
paths are provided. 

Random walks and tridiagonal matrices 

There is a very close relation between random walks and tridiagonal matrices, 
and the latter are known to be a powerful tool for counting walks with various 
specifications. Tridiagonal matrices are frequently used to describe the motion 
of a particle in a one-dimensional chain. 

Since the walks we are considering are made of steps ±1, we actually consider 
bidiagonal matrices of size N x N 



M(a,b) = 



01 



V 



61 



a.2 



62 




\ 



&JV-1 





(1.1) 



The starting point of our discussion is the fact that the explicit expansion of 
M(a, b)?j immediately suggests to order the various terms in the sum according 
to some property of the random walk that connects site i to site j in p steps. 

In more detail, if we identify the matrix element M(a, b)y with the sym- 
bol and agree that repeated neighboring indices are summed, we have: 
(M p )ij = (i, h)(ii, 12) ■ ■ ■ (i p -i,j). Since the matrix is bidiagonal, either ik+i = 
ik — 1 or ik+i = ife + 1. Because of this, the pair (ik, ifc+i) can also be identified 
with a move from site ik to a neighboring site ik+\ in a simple walk. Therefore, 
each sequence of p factors .. . (i p -i,j) not only represents a product of 
matrix elements, but also a random walk of p steps that connects site i to site 
j. The corresponding product of matrix elements may be interpreted as the 
"weight" of the walk. 

In this picture, the evaluation of M(a, b)^ consists in summing the weights 
of all walks of p steps from i to j. If we sum over lengths we get the resolvent, 
which is the generating function for all weighted simple paths from i to j: 

F(z, a, b; i, j) = zPM ^ h )h = t J - * M ( a > b )] Tj = 
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( dct[J — .zAf (a, b); j, z] 

1 ' det[J- zM(a,b)] 



(1.2) 



Here [A; i, j] is the matrix A with row i and column j removed. This important 
result is very well known (see for instance theor.4.7.2 in 

It is very useful to note that the matrix M(a, b) is similar to the matrix 
M(l,x) where X{ = dibi, 



M(l,x) = 



1 



V 





1 



•J'2 





1 



XjV-l 





(1.3) 



In the relation M(a, b) = SM(l,x) 5 1 the matrix 5 is diagonal, with entries 
s% = 1, S2 = ai, ... , sjv = «ia2 • ■ ■ ajv-i- Therefore 



[M(a,b)*'] ii = -i[M(l,xy] <3 

O 7" 



(1.4) 



From this point, we exploit the similarity ( |l.4| ), where we note that the 
factor Si/sj does not depend on p. To consider the matrix M(l,x), as we 
shall do hereafter, is a great simplification. In each walk, only upward steps 
(fc, k + 1) correspond to nontrivial factors Xk- Several walks may have the same 
weight since they contain, but with different order, the same intermediate steps 
(k, k + 1). By collecting these equal contributions, we have a useful expression 
in terms of certain counting numbers of random walks. The simplest example 
is 



N-l 

[M(l,x) 2 *V= £ ^ N( ni ,...,n t )x?x?...x? (1.5) 

t=\ K,(p,t) 

Here N(nx, . . . ,n t ) is the number of Dyck paths that make rij upward steps 
(j — from site j — 1 to site j, j = 1 . . . t. Their length is 2p = 2{n\ + . . - + n t ) 
and t is the height, which is bounded by the size of the matrix. 
The multiple sums in n\, . . . ,n t > 1 such that n\ + . . . + rij = p are summarized 
as the sum on the compositions of p into t integers. We denote the set of 
compositions of p into t integers as n(p; t). 

This formula gives the explicit representation of M(l,x) 1 ^ as a polynomial in 
the entries of the matrix. 

In the end of sect. 3 we describe the general case [M(l,x) p ]ij which corresponds 
to the counting of random paths of p steps from site i to site j, without positivity 
restriction. In the case x = 1 each walk has unit weight and one obtains the 



5 



counting number of simple walks from site i to site j of length p, restricted 
in the lattice of points from 1 to N, in terms of a simple matrix quantity. 
Since eigenvalues and eigenvectors are known, the multiple sums of the counting 



numbers in eq.(3.49) are equal to the trigonometric sum: 



2 p+1 A „ / kir \ . /. kn \ . /. kn 



[M(l, l) p ]i j = > cos p sin i sin j— (1.6) 

1 y ' ' 1 J N+l^ \N+1J \N+lJ \ J N + 1J v ; 

One can use random bidiagonal matrices to count random walks with re- 
strictions. The trick is to assign the variables Xi a probability density such that 
the contribution of unwanted walks is made to vanish in the average. 
For example, in |Tl| we studied the "q-root of unity" matrix ensemble, gener- 
alizing a problem first suggested by Zee [pi. A random matrix M(l,x) in the 
ensemble is characterized by independent random complex entries Xk chosen in 
the set of the q — roots of unity, with uniform probability. All moments 
have values 1 or 0, according to r being a multiple of q or not. Therefore, only 
walks that visit each site a number multiple of q do contribute, and with weight 
one, to tr M(l,x) p . This non local constraint strongly modifies the statistical 
properties of the random walks, as it was shown Noh et al. || in a model of 
surface growth. The case q = 2 of "even-visiting walks" has been investigated 
in great detail by Bauer et al. [|o), while the limit case q = oo was considered 
by Derrida et al. ||. 

On the other hand, given an ensemble of bidiagonal random matrices M(a, b), 
the enumeration methods may be useful for the study of the spectral density, 
usually in the limit N — > oo. If all random numbers Xj = ajbj are positive, 
the matrices are similar to real symmetric ones, M(c, c) with Cj = y/xj. The 
spectral density has support on the real axis and can be achieved by the usual 
resolvent approach. The average resolvent can be obtained by evaluating the 
ensemble average trM(l,x) p for every p. To this end, the random walk de- 
scription could prove to be a valuable tool, even for approximate results, if one 
were able to acknowledge and count the dominant walks. In the case where 
the random variables unrestricted, as in ref. |ll| , the spectral density has 

complex support, and the evaluation of the resolvent may only allow the deter- 
mination of the boundary of the support. 



2 Counting by iterations along the length of the 
path. 

One is often interested in the enumeration of paths according to a parameter 
having the additive property: the parameter of a sequence of Dyck paths is the 
sum, or a linear combination, of the parameters of the single Dyck paths in the 
sequence, see for instance ref. [|l2| . Perhaps the best known example of additive 
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parameter is the number of returns to the origin, see ref.|26|]. The method of 
enumeration by iteration along the length is suitable for the case of additive 
parameters, and it is here described following the pattern of paper pffl. 



Statistics of peaks and valleys. 

Given a walk we say that an 'inversion' occurs at time k if /ik 7^ Mfc+i- One 
may further partition inversions into peaks (fik = 1 , Mfe+i — — 1) an d valleys 
(m* = — 1 i Mfc+i = 1 )• We now show that the evaluation of the number of 
inversions in simple random walks is a straightforward combinatorial exercise 
rather similar to the evaluation of the number of visits to the origin, which is 
known in the literature. The counting of inversions in long walks will be ex- 
pressed in terms of the counting in shorter and simpler ones. 

It is useful to code a simple random walk of n steps as a sequence of to 
positive integers ni...n m , that count consecutive equal steps. A walk with 
initial step fii = 1 and described by the sequence {n\, . . . , n m ) begins with n\ 
steps /i = +1, followed by steps in opposite direction /j, = — 1, and so on. 
The sum n\ + . . . + n m is the length n of the walk, and the number of inversions 
is m — 1. For instance, the walk of 12 steps and 3 inversions in Fig.l is coded by 
the sequence (3,4,2,3). The same sequence also describes a second walk, with 
opposite signs of steps. 

Each sequence (ni, ri2, ■■,n m ) of positive integers is a composition of the integer 
n in to parts and the number of such compositions is C(n, m) 



It follows that the number of random paths of n steps and to inversions is 
2C(n, to + 1). By summing over inversions one reproduces the total number of 
random paths with n steps, 2^ m C(n, to) = 2". 

Our first evaluation is the number of walks of given length and given num- 
ber of inversions. It is convenient to consider first walks returning to the origin 
(therefore the number of steps n is even) and introduce the counting numbers: 

- Ci(n,m) is the number of strictly positive closed walks of n steps and m in- 
versions, m — 1, 3, ..,n — 3 

- C2(n, m) is the number of weakly positive closed walks of n steps and m inver- 
sions, to = 1, 3, .., n—l 

- 03(71, to) is the number of closed walks of n steps and to inversions, with no 
positivity requirement, to = 1, 2, .., n — 1. 

Since for strictly positive walks it is fix = fj,2 = 1 and /i„_i = ji n = —1, and 
for weakly positive walks it iss [i\ — 1 and fi n = —1, it follows that 




(2.1) 
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ci(n,m) — C2(n — 2,m) , for n>2 ; 
ci(2,m) = ca(2,m) = 5 m ,i 



(2.2) 




Fig.2 

A walk contributing to Ci(n^ m) either is in the class contributing to Ci(n, m) 
or has at least one visit to the origin before the last site, (Fig.2). If a walk 
contributing to ci but not to c\ first returns to the origin after n\ steps, we 
have 



<"2 



(n, m) = c\ (n, m) + ci(rii,mi) C2(n — ni,m — mi — 1) (2.3) 



1=2,4,. .,n—2 
1=3,5, ..,Tii -3 



In terms of generating functions 

Ci( x ,y)= X! c i {n,m)x n y 7] 



(2.4) 



77=2,4,.. 
-7=1,2,. .,77-1 



the eqs.(2.2), ( |2 .3[ ) translate into 

Ci(x,y) = x 2 y + x 2 C 2 (x,y) ; 
C 2 (x,y) = Ci(x,y) + yCx(x,y)C 2 {x,y) 

The resulting algebraic quadratic equations lead to 



Ci(x,y) 



1 — x 2 + x 2 y 2 — — x 2 + x 2 y 2 ) 2 — Ax 2 y 2 
2y 



(2.5) 



(2.6) 



^ = l-x 2 - x 2 y 2 - - x 2 + x 2 y 2 ) 2 - Ax 2 y 2 
' 2x 2 y 

The generating functions are closely related to the function u{x, y) evaluated 
by Narayana pTfl 
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n — 1 \ / n — 1 



u(x.y) = 1 — yjl — x + a;y) 2 — 4a;y 

oo to— 1 2 

= i + iy + VV-xV 

'-^ n — ] 

n =2 r=l 

For example, it is simple to obtain 

f) = E E ^-v-"- 1 -^ f " " ( " " 1 

* — ' * — ' n — IV m I \n — m 

n—2 rn—1 v 7 v 



(2.8) 



(2.9) 



A walk contributing to 03(71, m) (an example is given in Fig. 3) either it is 
strictly positive or negative, or it has a first return to the origin after ri\ steps, 
where it may or may not have an inversion. 




Fig.3 

This leads to the equation 

C 3 (x,y) = 2Cx(x,y) + {\ + y)C 1 {x,y)C i {x,y) 



(2.10) 



The solution is found using the previous results, and is written in the convenient 
form 



C 3 (x,y) = -1 + 



1 - x 2 (l - y) 2 



y/(l - x 2 + x 2 y 2 ) 2 - Ax 4 y 2 



(2.11) 



To extract the coefficients of the power expansion we use the generating 
function of Legendre polynomials, also found in paper |27]] 



- x 2 - x 2 y 2 ) 2 - 4x 4 y 2 



5> 2I V-1)"P„ 



n>0 



E 



y 2 + 1 

y 2 -l 



x y 



(2.12) 



Then: 



C 3 (x,y) = 2x 2 y + 2 



„2» 



?;>2.m>l 



..2m— 1 



n — 1 



2m 



n-1 

m 



n — 1 
m — 1 



(2. 



Eg. ( 2. 15 ) has a simple interpretation when the paths are coded by the se- 
quences mentioned at the beginning of this section. A random path with n steps 
and to inversions contributing to 03(11, to) is coded by a sequence (n%, n%, .., n m +i) 



such that V. odd nj = £ 



j even J 



n/2. The two subsequences s dd — (ni, 713, ..) 



j j odd 

and s even = ("-2,^4, ■■) both have (m + l)/2 terms if to is odd, whereas if to 
is even s dd has y + 1 terms and s et , e „ has to/2 terms. Then if to is odd, the 
number of sequences (n\, 122, .., n m+ i) is given by C(#, ^ + l) 2 , whereas if to is 
even it is given by the product C(^,tt + 1) C(^, Since each sequence de- 



scribes two walks, one obtains the coefficients Cs{n,m) in eq.(2.13). Analogous 
results are quoted in the monograph p8|. 



Let us now consider random walks of n steps and to inversions without 
the restriction that the end point of the path is the origin, and no restriction 
about positivity. Such walks co rres pond to all sequences (ji\,ni, ..,n m+ i) with 
iii + • • ■ rim+i = n. From eq^^T) we obtain the generating function of the 
counting numbers: 



C 4 (x,y) = 2Y / x n y m ( n ^ = 7T- U1 2 7T— T] (2-14) 
^ \ m J (l + y)[l-x(l + y)] 

Finally we evaluate the number 05(71, to) of weakly positive walks of n steps 
and to inversions, which may not be closed; an example is given in Fig. 4. The 
corresponding generating function is C^,(x,y). 




We may determine 05(71, m) in terms of shorter walks, by partitioning the set 
of weakly positive walks in three disjoint sets : a) the set of walks that after an 



10 



initial upward step never return to level one, b) the set of walks with an initial 
upward step followed by a walk in the ensemble described by the generating 
function C2(x,y); c) the set of walks of type b followed by a generic walk of 
type Cs(x, y). This translates into the equation 



Cs(x,y) = x + xC 5 {x,y) + xC 2 (x,y) + xy C 2 (x,y) C 5 (x,y) (2.15) 
Then 



1 + x — xy 1 ^/(l — x 2 + x 2 y 2 ) 2 ~ Ax 2 y 2 

&5{x,y) = h = 

2 2 1 — x — xy 

= xy+ 1 + X 2 +Xy C 3 (x,y) (2.16) 

The last equality allows an easy derivation of the counting numbers cs(n, m). 
As a simple check, we notice that the function C${x,y — 1), enumerates walks 
of fixed length and any number of inversions: 



The numbers are well known, see for instance |26[ 



Statistics of the area. 



The approach illustrated in the previous paragraph can also be used to ob- 
tain equations for the generating functions that enumerate the positive random 
walks of T steps returning to the origin and enclosing with the time axis a fixed 
area A. This problem was recently discussed by Jonsson and Wheater, and we 
refer to their paper [^3j for the physics motivations of the study and the analysis 
of the functional equation for the generating function. 



Let c\{A, T) and c 2 {A, T) be the numbers of, respectively, strictly and weakly 
positive walks that return to the origin after T steps, such that the area between 
the walk and the hori zont al axis is A. In both cases T is an even integer. Just 
as in Fig. 2 and in eq.(2.3) we have 



c 2 (A,T) = Cl {A,T)+ c 1 (A u T 1 )c 2 (A-A 1 ,T-T 1 ) (2.18) 

Ai=l,2,-. 
Tl=2,4,.. 

Any weakly positive path may be elevated by adding a step at the beginning 
and another at the end, thus transforming it into a strictly positive path, then 
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c 2 (A,T) = ci{A + T+l,T + 2) for A>2 
c 2 (A=l,T) = Cl (A=l,T) = 5 T , 2 

We define the generating functions 

F l (x,y)=J2^(A,T)x A y T for i = 1 , 2. 



(2.19) 



(2.20) 



Then eqs.( |2.18| ), @.19| ) translate into 



F 2 (x,y) = F 1 (x 1 y) + F 1 {x,y)F 2 (x,y) ; 
Fi(x,y) = xy 2 + xy 2 F 2 (x,xy) 

One obtains the nonlinear functional equation of Jonsson and Wheater 

F\ {x, y) = xy 2 + F x (x, y) F x (x, xy) 

which has a formal solution as an infinite continued fraction 

2 



Fi(x,y) = 



xy 



1 



x 3 y 2 



1 _ T^T7 



In a similar way we obtain the functional equation for F 2 (x, y) 
l-xy 2 [l + F 2 (x,y)} ' 



(2.21) 



(2.22) 



(2.23) 



1 + F 2 (x,xy) 



(2.24) 



F 2 (x,y) = -1+1- 



xy 



(2.25) 



From cqs. ( ^22| ), ( pi| ), one has = 1 ^ 4y2 and F 2 (l, y) 



- 2y z 



As shown in |23J, eq.( 2.22 ) allows the evaluation of -§^F\{x, y)\ x =i, which is the 
generating function for the area associated to the set of strictly positive paths 
of length T 



d_ 

dx 



:1 = ^/^A C1 (A,T) = T 



V 



- V 



(2.26) 
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In the same way we obtain the generating function for the area associated to 
the set of weakly positive paths of length T 



Ox 



F 2 (x,y)\ a 



1 - 2y 2 - v / T^4p 
2y 2 (l - Ay 2 ) 



(2.27) 



The results ( 2.26|) , ( 2.27 ) were obtained in [Q and extended to higher moments 
of t he ar ea in p5[ by a different method. We do not proceed to the analysis of 
eq.(2.22), whose main properties were found in paper p3fl . In the next section 
we describe a different approach to enumerate random walks, which will provide 
a different viewpoint of this problem. 



3 Counting by insertions which increase the width 
of the path. 

The counting method described in sect. 2 was based on the construction of the 
ensemble of random walks by adding new random walks at the end of previously 
considered ones. We now present a different approach where the ensemble is 
generated by inserting new random walks on the vertices of previously consid- 
ered ones. 

To evaluate the number of Dyck paths of 2n steps, height t, that visit pre- 
scribed numbers Nj of times the sites j = 0, 1 ... t, it is useful to consider the 
one-to-one map between Dyck paths and planar rooted trees [Q. The corre- 
spondence is shown in Fig. 5 where the rooted planar tree on the right side is 
mapped to the compressed "mountain range" of the left side, where each edge 
of the tree arises from one matched up-down pair. 



j N, n, 

5 2 2 

4 4 2 

3 5 3 

2 5 2 

1 3 1 

2 



Fig.5 
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Given a Dyck path, let nj denote the number of up steps between level j — 1 
and level j. Since Nj is the sum of the number of up steps from level j — 1 and 
level j plus the number of down steps between level j and level j + 1, we have 

Nj = nj + rij+i , j = l,..,t ; n = , N = m + 1 (3.1) 



The following table recalls corresponding elements : 



Dyck path 


planar rooted tree 


one pair of matched up-down steps 
length of path=2n 
Nj , j = l,.,t 
nj , j = l,.,t 
Nj-nj+i , j = l,..,t 
height t 


one edge of the tree 
tree has n edges 

sum of degrees of vertices at level j 
number of edges between level j — 1 and level j 
Nj — nj+i — Vj — number of vertices at level j 
height t 



The number N(v\, v 2 , ■■, i>t) of planar rooted trees with given set of non-root 
vertices {v±, is evaluated in a recursive way because planar rooted trees 

with a set of non-root vertices {v±, .., v t ,Vt+i} are obtained by adding Vt+i new 
vertices at level t + 1 and using t + 1 new edges to connect all of them to some 
of the vt vertices at level t. This may be done in (^i^"* 1 ) ways. Then 



Since Vj = Nj — = nj the evaluation ( p.2[ ) provides the number 

N(ni, n 2 , -Tit) of Dyck paths of height t, with specified numbers nj of up steps 
between level j — 1 and level j 

ATI \ f nx+n 2 -l\ f n 2 + n 3 -l\ (n t -i+n t -l\ 

Such paths have total length 2n, where n = n\ + .. + n t . 

We remark that it might seem more natural to enumerate Dyck paths of 
fixed length in terms of the integers Nj enumerating the visits at each level, 
rather than in terms of the positive integers nj enumerating the up steps from 
each level. The second coding is superior because the positive numbers nj are 
independent variables, save the restriction Y^, n j — n - 



We obtained the result ( p. 3D by inserting new mountain ranges on top of the 
highest vertices of Dyck paths |ll[ . We were unaware of the previous works on 
the enumeration of vertices in planar rooted trees |l^] , J2l|] . 

We now proceed to discuss some properties of the enumeration numbers 



(3.2) by introducing generating functions . They have a very direct relation 
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with general tridiagonal matrices. This fruitful and interesting connection was 
described long ago by P.Flajolet, in his basic paper apparently almost ignored 
in later investigations p2j. 

Let F t {x\, X2, --Xt) be the generating function for the numbers N[n\, . . . , nt): 
F t (x-L,x 2 , ■ ■ • ,x t ) = ^2 N(m,n 2 , ■ ■ ■ ,n t )x^x^ . . .x™' (3.4) 

rij— l,.,oo , j=\,.,t 

We here derive various useful properties of the function, that originate from a 
recursive relation. 

Proposition 1 

F t (xi,X2, ...,x t ) = F t -\{xi,x2, • • • , ) - F t -i(xi,X2, ..,xt-i) (3.5) 

1 - x t 



Proof: 



2J N(ni,n2,...,n t )x" 1 

nj— l,.,oc , j — l,.,t 



n t + n t -i - 1 



X t 



^2 N(n 1 ,n 2 , ..,n t ^ 1 )x™ 1 x 2 l2 ..Xtl 1 1 ^ 

nj — l,.,oc , j — l,.,t — 1 n t — 1 

2 N(n U Ti2,..,n t -i)x^x^..x^- 1 1 [(1 - x t y n ^ - l] = 

rij— l,.,oc , j — 1,.,£ — 1 

Ft_i(xi,x 2 , .., 1 *~ X ) - i*t_i(xi, a; 2 , ..,»t-i) • 
1 - x t 



Beginning with A^(n) = 1, one iteratively constructs the first few generating 
functions: 

F 1 (x) = -^ ] F 2 (x,y)= 



Fz{x,y,z) 



1 — x '' (1 — x)(l — x — y) 

xyz 



(1 — x — y){\ — x — y — z + xz) ' 
They suggest the following formal solution to equation ( |3.5| ): 

Proposition 2 

XiX 2 ..X t _iX t 



F t (Xl,X2,...X t --L,X t ) - r-5 -. r- 

P t (Xi, ..,X t )P t _i(Xi, ..,X t -l) 

P t -l(x 2 , ..,X t ) P t -2(x 2 ,..,X t -l) 

(3.6) 



P t (xi,..,x t ) P t -l(xi, ..,X t -l) 
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where Pt(xi, ..,Xt) is the polynomial generated through the recursion 

P t (xi,..,x t ) = P t -i{xi, ..,x t -i) - x t Pt- 2 (X1, -,X t -2) 

P = l , Pi(x) = l-x (3.7) 

Proof: the first equality is proven by direct substitution of the solution into the 
recursive formula (3.5) and by using the following identity, which is proven by 
repeated use of eq. (|3.7|): 

P t (x 1 , ..,x t ) = (1 - x t )Pt-i(xi, ..,x t -2, t-1 ) (3.8) 

1 - x t 



To prove the second equality in eq.(3.6), we formally solve the recurrence relation 
for polynomials, with the given initial conditions, by means of a transfer matrix: 



Pt{xi,-,x t ) P t ^ 1 (x 2 ,..,x t ) 
P t _x(xi, ..,a; t _i) P t - 2 (x 2 , -,x t -i) 

1 -xi\(l 1 
1 J \1 

By taking the determinant of both sides, we obtain a useful identity: 






-x t \ 


I 


-x t -i \ 




o ) 




o ) 



(3.9) 



Pt{x\, ..,Xt)Pt-2{%2, -,X t -l) - P t -l(xi, ..,X t -l)Pt-l{x 2 , ..,X t ) = -X\X2 ■ ■ -X t 

(3.10) 

The equality follows by dividing both terms by Pt{x\, ..,Xt)Pt-i{x\, ..,x t -i). • 

Note that the recursive property implies that the polynomial P{ {x\, .., Xt) 
coincides with the determinant of a matrix of size t + 1 : 

P t (x 1 ,..,x t )=det[I + M(l,x)} (3.11) 

The generating function for walks with height not greater than t is 



$t(xi, ..,x t ) = l + F 1 (x 1 )+F 2 (x 1 ,X2) + ... + F t (x l ,..,xt) 
P t -i(x 2 , ..,x t ) 
Pt(xi, ..,x t ) 



(3.12) 



This function has a simple representation in terms of the bidiagonal matrix 
M(l,x), of size t + l: 

*t(xi,x 2 ,...,xt) = [/ + M(l,x)]^j (3.13) 



By replacing the last term in the sum (3.12) with the right side of eq. ( p.5[ ) and 
resumming, one obtains the recursive property 

®t(xi,x 2 , ■ ■ .,xt-i,x t ) = $ t -i{xi,x 2 , • • ■ , -{———) (3.14) 

1 - x t 



1G 



and, since <i>i(a;) = (1 — x) , the following finite continued fraction represen- 
tation: 



®t{xi,x 2 ,...,x t ) = Wl (3.15) 



1 - 



■'2 



1 - x t 



The numbers N(n\, ...,nt) bring a very detailed information about random 
walks. One is often interested in reduced information, which require restricted 
summations over the variables n\...nt. This is, in general, a rather formidable 
task. We here provide some interesting examples. In the first we count closed 
weakly positive walks of fixed height and length and obtain identities involving 
sums of products of binomials. In the second example, we again discuss the 
problem of the area, by relating it to certain tridiagonal matrices. The third 
example addresses the issue of sums of weighted walks. 

Example 1: enumeration of walks with fixed length and height. 

Let C(n,t) be the number of weakly positive closed walks of length 2n and 
height t. Evidently, the parameter t cannot be larger than n. The sum over 
possible t values gives the well known number of weakly positive walks of length 
2n: 

£C(M> = ^( 2 :) (3.16) 
t=l 1 \ / 

The generating function for the numbers C{n,t) can be obtained by setting all 
arguments Xi equal to x in the function (|3.4|): 



F t (x) = J2 x n C(n,t) = 



x l _ P t -i(x) P t -i{x) 
x Pt{x)P t -i{x) P t (x) P t _i{x) 

C(n,t) = ^2 N{ ni ...n t ), {n>t) (3.17) 

n(n;t) 

where, for shortness, we adopt the notation for the sum over the t integers 
rii > 1 such that ri\ + ... + n t = n as a sum over the set n(n; t) of compositions 
of n into t integers rii > 1. The polynomials Pt(x) are easily evaluated from the 
recursion property 

P t (x) = P t -i{x) ~ xP t -i{x\ P = 1, Pl(x) = 1-x (3.18) 

and have a simple expression in terms of the roots of the equation z 1 — z + x = 0: 



t+2 t+2 1 

Pt(x) = Z -± ^ = -J= 

Zi - Zi V 1 - 4x 



1 + vi - 4i\ t+2 /l-vr - !^ 2 * 
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[(*+l)/2] , . 

K7T 



fc=l 



J] (l-4xcos 2 ^) (3.19) 



For the last equality we used the known eigenvalues of the matrix M(l, 1). 
Pt(x) is related to a Chebyshev polynomial of the second kind. 



The multiple sums over the compositions in eq.(3.17) are obtained by extracting 
the coefficient of x n from the generating function F t (x). This is simple for small 
t; for example: 



C(n,2) = Yl 

re(n;2) 

C(n,3) = 2 

K(n;3) 



n 2 + m — 1 

"■2 

n 2 + ni - 1 \ fn 3 + n 2 - 1 
n 2 / \ n 3 



->n— 1 



3 + V5\ / 2\ / 3- V5\ / 2 \ 

, , 7= I I 



We now consider the general case. 
Proposition 3 

C(n, t) = - J- V + - J- V 1^ (3.21) 

V ' y f + 2 ^ a " +1 i + 1 ^ b n+1 

3 3 3 3 

where n> t, {aj} and {bj} are, respectively, the zeros of Pt(x) and P t -\{x). 
Proof: we use the following identity for the inverse of a polynomial P(x) with 
simple roots Xj : 

to obtain an expression for the inverse of the polynomial Pt(x)Pt-i(x) in terms 
of the known roots of the factors: 

1 1 1 1 1 



P t (x)P t ^(x) ^ PKa^Pt^iaj) x-a 3 ^ P t (b 3 )P' t _ x (bj) x - b } 
We then use the following results, which will be proven 

Piia^Pt-^a-) = (t + 2)y^-, Pt{bj)PU(h) = -(« + 1)t3^T ( 3 - 22 ) 
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and obtain the power expansion of Ft(x) we searched for: 



P t (x)P t -i(x) 



-E 

4- 9 



1 - 4a.,- a; 1 



a' x — a, t + 

3 3 3 



1 \p 1 - ' 



1 - 4&J 



J X - b 3 



1 >— ^ 1 - 4a,- 1 ^ 1 - 46,- 
"t + 2^ a™ +1 + t + 1 ^ 6™ +1 



t + 2 



E* n 



To derive the first identity in eq.( p.22| ) one evaluates P'{x) in terms of the roots 
Z\ and z 2 of the equation z 2 — z + x — 0. For x — aj it is 

t+i , t+i 
P/(a J -) = -(t + 2) Zl + \ 

(zi - z 2 y 

The condition Pt(aj) — implies that z* +2 = Z2 +2 i then: 



21+2 J2t+2 



Pl(a j )P t . 1 (a j ) = -(t + 2)^ 



z-i — z 



[zi - z 2 } 3 



= -(*+2) 



[zi.z 2 ] 



This, together with z\ + z 2 = 1 and 2:12:2 = aj, allows the simplification that 
yields the result. The second identity is proven along the same line.* 



The structure of the coefficients C(n, t) as a two term difference suggests to 
introduce the numbers C(n, < t) of walks of length 2n and height not larger 
than t: 



C(n,<t) = C{n,l) + C{n,2) + ...C{n,t) 



-E 

4- 9 



1 - 4a. 



t + 2 T flj 



n+1 



( + 1 



t + 2 



£ 4 



3=1 



l cos 



2 -7 71 " 



t + 2 



sm 



2 -7 71 " 



2 E 



2/7 

n + jf(t + 2) 



t + 2 



[ t + 2 ] 

E 

L t + 2 



(3.23) 



2n + 2 
V?i + l + £(t + 2), 

(3.24) 



where we inserted the explicit form for aj from eq. (|3.19| ). The form of eq.( 3.23 ) 
may be recognized from the property C(n, < t) = [M (1, l) 2 "]n, where M(l, 1) 
has size t + 1. Of course C(n, t) = C (n, < t) - C(n, < t- 1). Note that C(n, < n) 
is precisely the number in eq. ( 3.1 6| ) . 

The generating function <&t(x) for the numbers is obtained by setting all argu- 
ments equal to x in the function (3.13). 



$t(x)=Y,C(n,<t)x n = - 



t+i 



sm 



2 J7T 

t+2 



n=0 



2^1 

3=1 



4a; cos 2 



(3.25) 
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Example 2: enumeration of walks of fixed length and area. 

The area of closed weakly positive walks with numbers of upward steps n\, . . . , rit 
is A = n\ + 3n2 + 5ri3 + . . . + (2t — l)nt, while the total length is T = 2n, 
n = ni + . . . + tit- With the position x% — xy, X2 = xy 3 , Xt = xy 21 ^ 1 in 
equation (jO]) we obtain the generating function F t (x,y) of counting numbers 
C(n,t,A) of positive closed walks with fixed length 2n, height t and area A: 



F t (x,y) = 5>VC(M,^4) = 

n,A 



,xy 



Pt-2{xy 3 



,xy 



2t-3 



■ ,2t-3\ 



C(n,t,A) 



P t {xy,..,xy 2t l ) P t -i{xy,..,xy 
,n 2 ,...,n t ) 

fc(rx;t) k 



(3.26) 



where the S function imposes the restriction that the integers n\ . . . n t of the 
compositions of n should correspond to walks with fixed area A. 
If we relax the restriction on height, being interested on counting numbers 
C(n, A) for the length and the area, the generating function is 



$oo(x, y) = l + Y,C(n,A)x n y A = l + J2Ft(x,y) 

n,A t>l 

Poo(xy 3 ,xy 5 , . . .) 
Poc(xy,xy 3 , . . .) 



(3.27) 



We shall solve the counting problem by giving first the explicit expression 
of the polynomial Pt(x,y) = Pt(xy 1 xy 3 , . . . ,xy 2t_1 ), which coincides with the 
determinant of the tridiagonal matrix of size t + 1: 



I + M(l,xy) = 



/I xy 
1 1 
1 



1 



xy 



1 1 xy 1 



(3.28) 



Proposition 4 

Pt(x,y) = 



[(*+l)/2] 

J2 i.~l) n x n y n{2n - 1] Yl 

n=l k=l 



1 



y 



i / 



2(t-n+2-fe) 



1 



,,2k 



(3.29) 



Proof: the polynomials solve the recurrence relation P t (x,y) = P t ~i(x,y) — 
X y 2t ~ 1 P t _ 2 {x,y) with Pq — 1 and P\(x,y) = 1 — xy. The explicit computa- 
tion of the first few ones suggests the following general structure: Pt(x,y) = 
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Y^n=a^ 2 ^ ( — x ) nc n,t{y), where c n j(y) are polynomials for which recurrence re- 
lations are easily written. By solving the first few ones, one easily arrives to 
conjecture the form stated in the proposition. A formal proof then follows by 
induction, but we omit it.». 

For the purpose of the statistics of the area, we need the limit t — *■ oo: 

n(2n-l) 

P 00 (x,y) = l + ^(- a; ) n - y - — (3.30) 

^ (i-r)(i-y 4 )---(±-r™) 

From this expression, we derive the following functional equation: 

Poc(x,y) - Poc(xy 2 ,y) + xy Poo(xy 4 ,y) = (3.31) 

Proof: 

n(2n— 1) 2n 

^ (i - y ) ■ ■ ■ (i - v ) i - v 

= Pooix, y) - }j- x r (1 _ y2) (1 _ y2( „_ x)) = 

= Poo(x,y) + xy P OQ {xy A ,y) • 

Since $00(2;, J/) = Pooixy 2 , y)/P OQ {x, y), we obtain a functional equation which 
corresponds to eq. (2.24) for F 2 (a;,?/): 

$oo(.T,y) = 1 + xy^> 00 {x,y)^> 00 (xy 2 ,y) (3.32) 



The formula 3.27 allows to compare $oo (a;, y) with the generating function 
F2{x 1 y) discussed in sect. 2, for weakly positive walks. There, the variable y 
was conjugated to the total length T = 2n and x was conjugated to the area. 
Here, for uniformity with the rest of the paper, x is conjugated to n and y to 
A. It turns out that ^(x 2 , y) = 1 + F 2 {y,x). 

Re mark. From the continued fraction representation of $00(2:1, x 2 , 2:3, ...), 
see eq.( 3.15 ) 

$ 00(2:1, 2:2, £3...) = ^ (3.33) 



1 Z2 



1 - 



1 - . . . 

it is easy to see that it solves the formal equation 

$00(2:1, x 2 , x 3 , .. .) = l + xx $00(2:1, 2:2,2:3,...) $00 (2:2, 2:3,2:4,...) (3.34) 

Let us suppose that each Xk is a function of two variables Xk = Xk{x,y). 
The above formal equation will turn into a functional equation for a function 
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$00 (x, y) if Xk(x, y) = Xk-i(x, y), where the replacements x — > x and y — > y do 
not depend on This is possible for the case 

Xk — x y pk+r , p and r real (3.35) 



Then eg. ( 3. 34 ) becomes 

$oo (z, y) = 1 + ^ (or, y) $oo(zy p , y) (3.36) 

Example 1 in this section corresponds to p — and r — , Example 2 cor- 



responds to p = 2, ?* = — 1. However the more general equation ( 3.3£ ) does 
not add substantially to the previous examples since it provides the generat- 
ing function for counting Dyck paths with given number of steps and given 
f(p,r) = ^2j(pj + r)rij, which is a linear combination of the area and the length 
of the paths: f(p, r) = {p/2)A + (2n)(r/2 + p/4). 
Another interesting case is 

x k = xy a " (3.37) 
and will be discussed in the end of the next example. 

As this work was being written we became aware of recent papers which have 
some overlap with this section j3(J J|l| . 

Example 3 : Toy functional integrals. 

In the introduction we mentioned the subject of weighted Dyck paths and spec- 
ulated about the possibility of an operative definition of path integral in the 
frame of Dyck paths. One may consider a variety of functionals G[y] where 7 
is a Dyck path of length 2n, and write 



C(n) 

where {sfc}|™ , with so = s 2n = 0, is the sequence of sites visited by the path 
7 and C(n) is the cardinality of the set T(n) of Dyck paths with 2n steps. It is 
given by the Catalan number C(n) — ^-j- ( 2 ™) , see eq.( |3.16 ). Since we shall 



not consider the continuum limit, our finite sums are here called toy functional 
integrals. 

We are also interested in the set T(n, t) of Dyck paths of length 2n, height t and 
cardinality C(n,t), see Proposition 3, whence one obtains the self-explanatory 
set r(n, < t). We discuss with some detail the class of functionals G[y] = 
J dry(7(r)) and consider the toy integral 

In 

Jl {n ^> J y ' ' 7 er(n,t) k=i 
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Two interesting examples are 

g (s) = s r , j(s)=c« s , a>0 (3.39) 

In particular, the toy integral with g(s) — s evaluates the average area enclosed 
by Dyck paths of height t and length 2n. The second example leads to explicit 
equatio ns, a nd will be discussed after the general case. 

Eq.( 3.3§| ) is first written in terms of the numbers Nj of visits of site j of the 
Dyck path 7, next in terms of the familiar numbers of upward steps rij 

In t t t 

Y 9( s k) = Y 9(J) N J = Y 9(M n J + n J+i) = Y^U) + 9(j - 1)K 

k=i 3=0 3=0 3=0 

(3.40) 

The sum over paths is performed by introducing the counting numbers N(rii, n 2l .., 
t t 

Y Y^(j)+9U-l)]nj = Y N{n 1 ,n 2 ,..,n t )Y[9(j)+9(j- 1 )}n j = 

jer(n,t)j = l K(n;t) j=l 

= ^G-C(n,t,G) (3.41) 

Gee 

where C(n, t, G) is the number of Dyck paths of length 2n, height t and fixed 
value G 

C(n, t,G)= Y N( ni ,n 2 , ..,n t ) 6 (g - + g(i - l)]n^ (3.42) 

and the sum G £ Q runs over the set Q of the possible values of G. The 
generating function 

F t (x,y) = YC(n,t, G)x n y G 

n,G 



is obtained from the general theory by setting Xk = x y9( k )+a(k-i) m e q.(3.6). 
It may be evaluated from ratios of the determinants Pt(x,y). 
The function 

$oo(x,y) = l + F 1 (x,y) + ... + F t (x,y) + ... = Y C ( n > G)x n y G (3.43) 

n,G 

is the generating function for the numbers C(n, G) = J2t=i t, G) that count 
Dyck paths of length 2n and fixed value G. The generating function f(x) for 
the numbers G n = ^ G< zgGC(n,G) can be written as the first derivative of 



d 

f(x) = Y GnxU = -ff®oc{x 7 y) 



(3.44) 
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While we have been unable to find a "quadratic" functional equation for the 



function &oo{x,y) related to the toy integral (3.38) with g(s) — s r for r > 2, 
we were successfull with the exponential function g(s) = ca s . In this case the 
functional evaluated on a Dyck path is, see eq. (3.40), 



fc=0 



a + 1 x - 

c ^a s 



a 



s=l 



The generating function for counting numbers is obtained with the position 
Xk = xy( a \ so that 

~ni rr, n t _ rr ,n 1 +...n t ma+Tl2a 2 +—+nta t _ n G 



V 



x y 



and, for the general theory, it is given by the continued fraction 

1 



*oo(x,y) 



1 - 



xy 



xy 



1 - 



xy 



i 



which corresponds to the functional equation 

$oc(x,y) = 1 + xy a <P 00 {x,y)<P 00 (x,y a ) 



(3.45) 



From this equation, since $00 (^, 1) = (1 — \/l — Ax)/( 2x), one easily obtains the 
generating function for the numbers G n , see eq.( 3.44[ ), 



d_ 

dy 



$oo(x,y) 



(1 - 2x - VI - 4a;) 



x (1 - a + (1 + a) VI - Ax) 



It is interesting that the counting numbers C(n,A) , see eq.(3.27), allow to 
define a partition function for Dyck paths of any length L = 2n, weighted with 
their length L and area A: 

[ £> 7 e-^m-bAh] = J2 e - aL - hA C(L/2, A) = $oc(e" 2a , e _b ) (3.46) 
Jr LA 

where T is the set of all Dick paths and a different normalization has been used. 

Enumeration of closed paths. 

The counting numbers N(ni, . . . , n t ) apply to weakly positive closed paths, and 
are useful for expressing the matrix clement (1,1) of Af(l,x) 2p as a polynomial 
in the entries Xi of the matrix, see eq. (II. 
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The generalization of the counting numbers to arbitrary closed paths is straight- 
forward, and is given in 
Proposition 5 

The number of closed paths of height t, depth s and length 2p, containing n_j 
upward steps (— i — i = s — 1, ... ,0, and nj upward steps (j — 

j = 1, . . . ,t, is 

N s , t (2p; n- a+1 ...n \n 1 ...n t ) = 
= N(n , n_i . . . ( n ° + ni \ N(n 1 ,n 2 , . . . ,n t ) = 

V m ) 

n 1 / n-i + n-i+i - 1\ / n + "A yr ( n t + n t+l - 1\ 
»-« A -o n l+1 ) (3 ' 47) 

Proof: Closed paths contributing to the above counting number are sequences of 
m closed positive paths and n closed negative paths. The factors N(n 0} . . . , 
and N(n±, . . . , n t ) respectively count the negative and positive closed paths that 
are obtained by joining the negative or positive subsequences. The intermediate 
binomial factor counts the ways into which the subsequences can be ordered. • 
The length 2p is specified in the symbol in view of a further generalization to 
open paths; for closed paths p = n_ s+ i + . . . + n t . 

The counting numbers allow to write the expression for the matrix element (i, i), 
i = 1 ... N of the powers of the bidiagonal matrix M(l,x) as a polynomial in 
the entries: 

N-i N-l 

Af(l,x)?? = 2 2 JV s ,t(2p;n_ a+ i...n |ni...n t ) 

i=0 s+t=l K (2p,s+t) 

x^r-- (3-48) 



Enumeration of general paths 

The enumeration is extended to paths that may not return to the origin. For 
open paths the counting numbers N Sj t(p ; n_ s +i, . . . no | ni, . . . n t ) correspond to 
an expression which is more general than the above one. 

The paths have length p and make the specified set of upward steps between 
sites — s and t. Since the total number of upward steps is u = n_ s+ i + . . . + n t 
and the number of downward steps is d = p — u, the endpoint of the walks is 
site (level) y = u — d = 2u — p. If 2m — p^O the paths are open. 
For all paths it is m,..,n t > 1, because the walk has to reach the highest site 
t. However, for open paths some numbers nj with i < can be zero. More 
precisely, if y > all m > 1 because the walk has to ascend from — s to the level 
y; if y < it may be that nt = 0, for k < y. 
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With this notation and postponing the problem of the evaluation of the 
counting numbers, we write a general expression for the power of a matrix 
M(l,x) as a polynomial in the entries x\, . . . , xn-i- 

Proposition 6 

[M(l, xfk, = £ J2 N '*to n -*+i ■■■no\ni... n*)^ 1 • • ■ x^x? ■ ■ ■ ^ t -i 

s,t n_ s -f i...rc.£ 

(3.49) 

In the first sum s and t have the following restrictions: 

1) the finite size of the matrix, which constrains walks inside the strip of integers 
1 . . . N — 1 that label the entries. For this reason: < s < i — 1 and < t < 
N - 

2) walks have to visit all sites between i and j included. Therefore, if j > i it is 
t > j — i; if j < i it is s > i — j. If i = j both s and t cannot be zero. 

3) the sum on n_ s+ i . . . n t has the restriction that, being u — m s + • ■ ■ n t the 
number of upward moves, the constraint 2u = j — i + p must be satisfied. (Note 
that for a bidiagonal matrix [M p ].;.j = if p and \j — i\ have opposite parities). 

The general expression for the counting numbers has been given by Krat- 
tenthaler |Is|| , including the case of the null step. We here discuss the case of 
paths that, starting from the origin, only visit sites i > 0. 

Proposition 7 

The number of positive walks of length p, that never return to the origin and 
with prescribed numbers of upward steps is: 

N , t (p; l,n 2 ,...n t ) = l[( nj+nj - 1 ~' 2 )N(ri v ,...,n t ) (3.50) 

j=2 V n 3 1 J 

where y = 2u — p is the height of the walk, u = 1 + 712 + • • • nt is the total 
number of upward steps (see comment in proposition 5). 

To obtain this counting number it is convenient to generalize the planar 
rooted tree by using edges of two types, with solid or dotted lines, as in Fig. 6. 
There is a one-to-one map between the "open mountain range" and the cor- 
responding generalized planar rooted tree: to any solid line edge of the tree 
correspond, as before, a pair of matched up-and-down sides of the mountain 
range, while a dotted line of the tree corresponds to a up side of the mountain 
range without the matching down side. 
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Fig.6 

The correspondence between the number of visits at each level and the number 
of up steps is 

Nj = rij + n j+1 - 1 , 1 < j < y 

Nj = rij + nj + i , y<j<t (3.51) 

Each vertex of the planar rooted tree is connected to a vertex at the lower level 
by a full line or a dotted line, rij counts the full and dotted lines in the tree 
between level j — 1 and level j; then 

Vj = n. } , j = l, ..,t (3.52) 

The number of planar rooted trees is evaluated first by drawing the root 
and a dotted path from the root to level y, with one vertex at each level 1, .., y, 
next by adding a number Vj of new vertices at each level j, and connecting 
each vertex with a full line to a vertex at the lower level. If we call Vj — Vj for 
j = y + l)..,t and Vj = Vj — 1 for j = l,..,y, we obtain that the number of 
planar rooted trees with given set of non-root vertices: 



N(v 1 ,v 2 , ..,v t ) 



v t + v t -x 



V t -1 + W t _2 
V t -1 



V 2 + Vi 
V2 



(3.53) 



From this, eq.(3.52) and the relation between Vj and Vj we obtain the number 
in proposition 7. • 
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